Tunneling anisotropic magnetoresistance in single-molecule magnet junctions 
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We theoretically investigate quantum transport through single-molecule magnet (SMM) junctions 
with ferromagnetic and normal-metal leads in the sequential regime. The current obtained by means 
of the rate-equation gives rise to the tunneling anisotropic magnetoresistance (TAMR), which varies 
with the angle between the magnetization direction of ferromagnetic lead and the easy axis of SMM. 
The angular dependence of TAMR can serve as a probe to determine experimentally the easy axis 
of SMM. Moreover, it is demonstrated that both the magnitude and sign of TAMR are tunable by 
the bias voltage, suggesting a promising TAMR based spintronic molecule-device. 

PACS numbers: 75.50.Xx, 75.47.-m, 85.75.-d 
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In the past two decades, magnetoresistance in mag- 
netic tunnel-junctions has received much attention 
because of its strong dependence on the relative 
magnetization-orientations of two ferromagnetic (FM) 
electrodes, in terms of which the tunneling magnetoresis- 
tance (TMR) devices are developed for various applica- 
tions in magnetic sensors and information storage tech- 
nology P, Hjj. More recently the tunneling anisotropic 
magnetoresistance (TAMR) beyond the conven- 

tional TMR has been observed in the presence of spin- 
orbit interactions, which depends on the relative direc- 
tions between magnetization-orientation of the FM lead 
and the crystal axis. Particularly, the TAMR even ex- 
ists in tunnel junctions with a single magnetic electrode 
such as Fe/GaAs/Au ||, while the conventional TMR 
effect does not appear in this configuration. Therefore, 
the TAMR effect may lead to new spintronic devices with 
only one magnetic lead d-Q. 

Single-molecule magnets (SMMs) what we consider 
possess high spins (S > 1/2) and uniaxial magnetic 
anisotropy with an easy axis, which have potential ap- 
plications in molecular spintronics @. In recent experi- 
mental 0, H| and theoretical [t^-Hfl studies on electronic 
transport through magnetic molecules many fascinating 
properties have been found, such as complex tunneling 
spectra Q , negative differential conductance @, EH U | , 
Kondo effect [l2| , Berry phase blockade , and colossal 
spin fluctuations In particular, the spin-polarized 

transport through a SMM shows that the magnetiza- 
tion of SMM can be controlled by spin polarized cur- 
rent fiol. Il5 - 17 j. spin- bias [l8[ and thermal spin-transfer 
torque [19( . In addition, a pure spin-current generated by 
thermoelectric effects [20|, polarization reversal of spin- 
current 21 1, spin diode behavior [ll[ and spin filter effect 
22, 23| are also found. The conventional TMR effect in 
spin-dependent transport through a SMM in the sequen- 
tial, cotunneling and Kondo regimes is also investigated 
by Misiorny et al. 24) , where the magnetizations of FM 
leads are collinear with the magnetic easy axis of SMM. 
On the other hand, they also investigate theoretically the 
magnetic switching of a SMM coupled to two collinear 
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FIG. 1: (Color online) Schematic of SMM tunnel junction 
consisting of FM and normal metallic electrodes. (j> denotes 
the angle between magnetization of the FM electrode and the 
easy axis of the SMM. 



FM leads, in which the magnetic easy axis of SMM forms 
an angle with the FM electrodes 251 ] . 

In this paper, we study anisotropic effects of the tun- 
nel junction (see Fig. [lj, which consists of a SMM 
sandwiched between one FM electrode and one normal 
metallic (NM) electrode, where the FM electrode is not 
collinear with the easy axis of SMM different from the 
TMR case. This setup could be realized by the scan- 
ning tunneling microscope with a spin-polarized tip and 
a SMM on a nonmagnetic metallic substrate [ll|, 27 , 28[ . 
We are going to show that the transport current depends 
on the magnetization direction of FM electrode, resulting 
in the TAMR effect. 

The system is described by a Hamiltonian ToL 18[ . H = 
Hieads + H S mm + H T , in which the first term H teads = 
Eo=l,r Ek,r=± £akcl kT c akT is for non-interacting elec- 
trons in the magnetic leads with c^ k (c Q kr) being the 
creation (annihilation) operator of an electron in the 
lead-a (with a — L, R denoting the left and right leads 
respectively) of energy £ a k and spin index r = ±. In ad- 
dition, the spin polarization of the FM lead-a is defined 



as P a = {p a+ - p a -)/(p a+ + p a -), with p a+ (-) denot- 
ing the state density of the majority (minority) electrons. 
The model Hamiltonian of SMM is 



Hsmm = ed i d x + Udld+dld- - Js ■ S - K (S z ) 2 



(1) 



2 



where d x (d x ) is the relevant electron creation (annihila- 
tion) operator in the lowest unoccupied molecular orbital 
(LUMO) level and s= J2 XX ' 4( CT xx' / 2 ) d x' is the 

corre- 
sponding electron spin operator (er is the vector of Pauli 
matrices, x — i denotes the spin index), e is the single- 
electron energy of the LUMO level, which is tunable by 
the gate voltage. U represents the Coulomb energy of 
two electrons with opposite spins in the LUMO level. J 
is the exchange coupling parameter between the spin-S 
of SMM and electron spin-s, which can be of either FM 
(J > 0) or antiferromagnetic (AFM) (J < 0) type, and 
the parameter K > describes the easy-axis anisotropy 
of SMM. Moreover the total spin is defined as S t = S + s, 
and the many-body states of SMM and electron can be 
written as \n, S t ; m), with n denoting the charge state of 
the molecule, St the total spin quantum- number, and m 
the eigenvalues of Sf , while the corresponding eigenenegy 

i s £ \n,S t ;m)- 

The tunneling processes between the molecule and 
leads can be described by the Hamiltonian [T|| TH 25 1 

H T = Eafcxr^fcrC 008 ^^ + XSin^^)^ + H.C., 

where t a denotes the tunnel matrix element of the 
molecule and the lead-a, and <p a (with the angle defini- 
tion = <fi and <pn = 0) is the angle between the magne- 
tization direction of lead-a and the easy axis of the SMM 
(as z-axis) . The spin-dependent tunnel coupling-strength 
is described by T a ± — F Q (1 ± P a )/2 for spin-majority 
(upper sign) and spin-minority (lower sign) electrons in 
the lead-a, with r a = F Q+ + r a _. 

For weak coupling between the SMM and leads, i.e. 
r Q <C fcfiT and |J|, the molecule has enough time to 
relax to the eigenstates of Hsmm between two consec- 



utive electron tunneling processes 18j. Therefore, we 
adopt the rate equation approach to study the station- 
ary transport in the sequential regime. Denoting Pi (t) as 
the occupation probability of the SMM in the molecular 
state \i) at time t, we have 



dt 



Y, w « ip « -w^' Pi- 



rn 



Applying the Fermi's golden rule T(J Hj], the sequential 
transition rate W^ 1 from the state \i) to \i'} with respect 
to the lead-a is given bv Til Til 25j 



w >/ = j2 ra{1 + h xPa) [\(i\b x \i>)f 

X 

f{e v -e i - l i a )^\b x \i)f f{t 



Si + fJ-a ) 



(3) 



where b x = cos 4f- d x 



X sin (i c 



= {-lf^eV/2, 
Ei being the energy of state \i) and f(x) is the Fermi 
distribution function. After some algebra, the transition 



rate Eq. Q can be rearranged as 

ur i,i' \ - T a (l+ X P a coscj) a ) , 2 

W c! =2^ ^p, l\(A d x\ l )\ 



f(si< -£i~ (ia) + \ d x \i)\ f{s v - e l + (i a ) 



(4) 



Furthermore, by introducing the effective tunnel strength 
T a ± = r Q (l ± P a cos (j> n )/2 with the effective spin po- 
larization P a cos 4> a [29( , the magnetization of FM lead 
may be considered as collinear with the easy axis of the 
SMM. 

The stationary probabilities obtained from the condi- 
tion, = 0, with the transition rate together result in 
the current through lead-a 

I a = (-l) 5 ^eJ2(^ - n v )W 1 / Pi, 

ii' 

and the TAMR ratio in the FM/SMM/NM tunnel junc- 
tions is defined as 151 



Rtamr{4>) 



no) - m 



(5) 



where I(<j>) is the current with the angle <fi between the 
magnetization direction of FM lead and the easy axis of 
the SMM. 

In Fig. Ufa), we show the bias- voltage dependence of 
the current for different angles <fi with the FM exchange 
coupling (J > 0), where only the angle range cf) £ [0, n/2] 
is considered, since for <j> e [tt/2,tt] one just reverses the 
polarization direction of FM electrode. The current spec- 
trum is asymmetric under the bias voltage reversal due 
to the asymmetric coupling of the SMM and the electron- 
spin in the FM lead. The current magnitude especially 
in the bias voltage region, V = 1 — 5 mV is more sensitive 
(with large changes) to the angle variation and increases 
with it. The current plateau decreases with the angle 
increase around V — —2 mV, while the situation is just 
opposite around the bias voltage value V = 2 mV. It is 
a well known difficulty that the easy axis of SMM rela- 
tive to the laboratory frame can not be well determined 
[26| in SMM transistor experiments 0, IH , and thus the 
angular-dependent current with respect to the FM lead 
makes it possible to detect experimentally the magnetic 
easy axis of SMM jH [n|]. 

The TAMR ratio as a function of the angle cj> for differ- 
ent bias voltages is plotted in Fig. Efb), which exhibits 
a periodic behavior with increasing the angle (f>. The 
TAMR is positive only for the bias voltage V = —2 mV 
(solid line), where the current decreases monotonously 
with the increase of the angle (f> [see Fig. [U(a)], and is 
negative otherwise. The TAMR amplitudes with V = 2 
mV (dash line) and V — 4 mV (dotted line) are larger 
evidently. Furthermore, when the bias voltage is high 
enough all transport channels open up and the TAMR 
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FIG. 2: (Color online) (a) The absolute value of current as 
a function of the bias voltage V for different angles (/>, (b) 
TAMR as a function of the angle (f> for different bias voltages 
V with parameters S = 2, e = 0.9 meV, J = 0.4 meV, U = 1 
meV, K = 0.05 meV, k B T = 0.04 meV, P L = 0.8, Pr = 0, 
r = r L =r H = O.OOl meV and I max = eF/h « 0.25 nA. 



becomes symmetric with respect to the positive and neg- 
ative bias voltages V = ±6 mV. 

Figure[3]displays transport characteristics for the AFM 
(J < 0) exchange interaction, in which the SMM state 
|l,3/2;m) has a lower energy than the state |l,5/2;m). 
Different from the FM case, the current plateau around 

V = ±2 mV decreases (increases) with the increasing 
angle <j> [Fig. and the TAMR for V = ±2 mV is 
positive (negative) [dash (solid) line in Fig. EJb)]. The 
TAMR amplitude has the largest value at the bias voltage 

V = — 2 mV and exhibits the same behavior as the FM 
case for V = ±6 mV. 

In order to understand the characteristic current- 
spectrum and the TAMR effect, the main probability dis- 
tribution of the steady molecular states as a function of 
the bias voltage V is plotted in Fig.|U where the FM elec- 
trode is collinear with the easy axis of SMM (i.e. = 0). 
We can find that the steady transport for positive bias 
voltages is determined by the states with negative eigen- 
values of spin operator S*, while the situation is just 
opposite for negative bias voltages. This is because of 
the spin-flip process between the transport electron-spin 
and the SMM 0, El Q E3 ■ 

In the FM case, the main occupied states of the 
SMM are |0,2;-2) and 1, 5/2; -5/2) at the bias volt- 
age V — 2 mV [see Fig. H[a)], and thus the transport is 
dominated by the transition 0,2; -2) |1, 5/2; —5/2) 



FIG. 3: (Color online) (a) The absolute value of current as 
a function of the bias voltage V for different angles (/>, (b) 
TAMR as a function of the angle <f> for different bias voltages 
V. 



of a spin-down electron. The effective tunnel strength 
= Tl (1 — Pl cos <fi) jl increases with the angle <f> vari- 
ation from to 7r/2, and therefore the current increases, 
while the TAMR is negative [dash line in Fig. EJb)]. 
On the other hand, when V = —2 mV, the transport 
is dominated by the transition [0,2; 2) 1, 5/2; 5/2} of 
a spin-up electron. Since the effective tunnel strength 
Tl+ = Tl(1 + Pl cos c/>)/2 decreases with the increas- 
ing angle 4>, the TAMR is positive [solid line in Fig. 
[5£b)]. At the bias voltage V — 4 mV, the additional 
main transitions |0, 2; —2) <^> |1, 3/2; —3/2} (spin-up) and 

1,3/2; -3/2) |2,2;-2) (spin-down) take part in the 
transport processes, and the TAMR exhibits a similar 
behavior as the case of V — 2 mV. While at V = —4 mV 
the TAMR amplitude is smaller due to the competition 
between different transitions of opposite spins. 

The state occupations of the AFM exchange interac- 
tion are plotted in Fig. 0Jb), from which we see that 
the transport is dominated by the transition 0, 2; — 2) <^> 

1,3/2; -3/2) (|0,2;2) 1, 3/2; 3/2)) of a spin-up (spin- 
down) electron at bias voltage V = ±2 mV, and thus 
the TAMR is positive (negative) [dash (solid) line in 
Fig. Eljb)]. The competition between different transi- 
tions of opposite spins suppresses the TAMR amplitudes 
as shown by dotted and dash-and-dot lines respectively 
for V = ±4mV in Fig. Ifb). 

In summary, the TAMR is obtained explicitly from 
the angular-dependent transport in FM/SMM/NM mag- 
netic tunnel junctions by means of the rate-equation ap- 
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FIG. 4: (Color online) The main probability distribution of 
molecular eigenstates as a function of the bias voltage V for 
FM (a) and AFM (b) exchange couplings in the collinear con- 
figuration (4> = o). 

proach. It is demonstrated that the angle between the 
FM lead and the easy axis of SMM is crucial to induce the 
TAMR through the exchange coupling between the SMM 
of uniaxial magnetic anisotropy and transport electron- 
spin. Both the magnitude and sign of TAMR are effi- 
ciently controllable by the bias voltage, suggesting that 
this SMM tunnel-junction can be a promising spintronic 
device. On the other hand, this TAMR effect may serve 
as a probe to detect the uniaxial magnetic anisotropy of 
the SMM. 
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